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Abstract 

Considering a linear parabolic stochastic partial differential equation driven by impul- 
sive space time noise, 

dX t + AX t dt = Q 1/2 dZ t , X =x eH, ie[0,T], 

we approximate the distribution of Xt- (Zt)te[o,T] i s an impulsive cylindrical process and 
Q describes the spatial covariance structure of the noise; Tr(A~ a ) < oo for some a > and 
A^Q is bounded for some /3 e (a - 1, a]. 

A discretization (-X^) n€ {o,i,...,jv} is defined via the finite element method in space 
(parameter h > 0) and a 0-method in time (parameter At = T/N). For ip e C 2 (H;&) we 

^-E^AV 1 ' " 

E<p(x£) - E<p(X T )\ = 0(h 2 ~< + (At) 7 ) 



show an integral representation for the error \Eip(X^) - E<,c(At)| and prove that 



where 7 < 1 - a + /3. 
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1 Introduction 

In this paper, we study the weak order of convergence of numerical approximations of the 
solutions of a certain class of linear parabolic stochastic partial differential equations (SPDEs, 
for short) driven by impulsive space time noise. Unlike the strong order of convergence which 
measures the pathwise approximation of the true solution by a numerical one (cf., e.g. 
|12j , |13] , |34] , [33] ) , the weak order is concerned with the approximation of the law of the true 
solution at a fixed time. There are not many works in literature about the weak approximation 
of the solutions of SPDEs (see [2], [8], [9], |10| . |14j ) and, to our knowledge, only SPDEs driven 
by Gaussian noise have been considered in this context so far. This work extends the paper 
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[9j by A. Debussche and J. Printems, where the following Hilbert space valued stochastic 
differential equation is considered: 

dX t + AX t dt = Q 1/2 dW t , X = x £H, te[0,T]. (1) 

Here A : D(A) c H -> H is a unbounded strictly positive definite self-adjoint operator whose 
domain D(A) is compactly embedded in H; Q : i7 -> i/ is a bounded nonnegative definite 
symmetric operator and (W^^mr] is a cylindrical Wiener process on H, T e (0,oo). A 
standard reference for this setting is |6]. 

If we set H := L 2 (C) = L 2 (O,B(O),d0, O c H d open and bounded, and := 
(-A, H 2 (0) n Hq(O)), then ([I]) is an abstract formulation of the stochastic heat equation 
with Dirichlet boundary conditions 



- AX(t,0 = v(t,0, (U) 6 [0,T] x O, 

X( v ) =0 on [0,T] xdO, 

X(0,-) = x on C 



(2) 



Here r) = |2 is a generalized random function which can be described as the time derivative 
(in a distributional sense, cf. |29], |33| ) of a real-valued generalized Gaussian process, formally 
written as 

V(t,0 = f 1o(t;,OW(t,Od(, (3) 

where (W(t, •))t€[o,T] i s a cylindrical Wiener process on L 2 (0) and go is a (generalized) 
function on O x O. The operator Q is given by Qx(£) = / g(£, 0^(0 <^C with g(£, £) = 
Jo ?o(C) r )<7o( r > C) ^ r describing the spatial correlation of the noise 77, cf. [27], Ch. 4.9.2 and 
Example 14.26. 

Throughout this article, let if := L 2 (0). We consider the equation 

dX t + AX t dt = Q l/2 dZ t , X = x £H, t e [0, T], (4) 

where A and Q are as above and (-^t)t £ [o,T] = (^(^> ' ))te[o,T] i s an impulsive cylindrical process 
on If, see Section 2 for the definition. This is an abstract version of problem (|2]) if one replaces 
W(t, C) by Z(t, in the formal definition (|3]) of the noise 77 and if furthermore go is symmetric 
and positive semidefinite in the sense that j f qo({,, C)V ; (OV ; (C) d£ dC, > for all test functions 

In [JS], a discretization (^^) n 6{l,...,JV} °f the solution (Ai)f £ [o,r] of equation ([1]) is obtained 
by the finite element method in space (parameter h > 0) and a ^-method in time (parameter 
At = T/N). Under the assumption that A~ a is a finite trace operator for some a > and that 
A^Q is bounded for some j3 e (a - 1, a], it is shown there that for functions e C 2 (H, R), 

lE^Xf ) - E^(X r )| < C ■ + (At)T) (5) 

for any y<l-a + /3<l. 

In this paper, we consider the analoguous discretization of the solution of (|4]) and make 
the same assumptions on the operators A and Q as in |9]j. We give a representation formula 
for the error (Theorem [5]) and, under some integrability condition on the jump size intensity 
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v of the cylindrical impulsive process {Zt)t<=[o,T]) we show that ([5]) holds also for the solution 
(Xt)te[o,T] m @ an d the corresponding discretization (Theorem |SJ) . 

SPDEs driven by impulsive noise (or Poisson noise) have been considered, e.g. in [I], |15| . 
|f 9j . |22j . |25j, [26J. The monograph |27| gives a good overview about SPDEs driven by Levy 
noise. In |16| and |17] . numerical approximations in time and space of SPDEs driven by Poisson 
random measures are investigated and the strong error is estimated. Of course, this especially 
implies an estimate for the weak approximation error. A difference to our result is that we look 
at impulsive noise which is white in time and coloured in space whereas in [16] and [IT] a class 
of SPDEs driven by Poisson random measures which correspond to impulsive space time white 
noise is considered. Our motivation for this paper was to show that the techniques applied in 
[9] with respect to the cylindrical Wiener process also work for certain jump processes. 

The main technical difference between fl4|) and ([1]) lies in the fact that the impulsive 
cylindrical process (Zt) t€ r 0T i is a purely discontinuous Hilbert space valued martingale, while 
the cylindrical Wiener process (Wt)te[o,T] i s continuous. As a consequence, the main tools for 
estimating the weak order of convergence for the numerical scheme — the Ito formula and 
(connected with it) the backward Kolmogorov equations for certain processes associated with 
the solutions of the SPDEs and their discretizations — are completely different for ((4]) and 
(H]). The main task therefore is to find manageable expression for the approximation error, 
which allows estimates using techniques similar to those in [9]]. 

We note that (|5]) remains true for the solution (X t ) t ^ 0T ^ of 

dX t + AX t dt = Ql /2 dW t + Q{ /2 dZ t , X = x tH, t e [0,T], 

and the corresponding discretization. Here, the covariance operators Qq and Q\ are assumed 
to have the same properties as Q above. 

Let us finally remark that the weak order for numerical approximations of ordinary stochas- 
tic differential equations with jumps has been studied in several papers, see [4j for a compre- 
hensive survey. In [28j the Euler scheme for Levy driven equations of the form dXt = f(Xf-)dZt 
is considered and it is proved that the weak order is O(At) if the jump intensity measure of 
Z has its first several moments finite. In [24] and |21] . the more general concept of stochastic 
Taylor expansions is used but the underlying jump intensity measure is supposed to be finite, 
i.e. only finitely many jumps appear on a bounded time interval. Comparing these finite di- 
mensional settings with (J4j) , we note the following: If H is finite dimensional the conditions 
on A and Q will be trivially fulfilled for a = f3 - 1 and one will get the same order as in |28] , 
Secondly, the jump intensity measure of the cylindrical process (^t)t e [o,T] with state space 
U d H - L 2 (0) is a measure on (U,B(U)) and has to be distinguished from the jump size 
intensity measure v as a measure on (R,i3(IR)), see Section 2. The conditions we impose on 
(•Zt)te[o,T] m order to obtain our results can be best expressed in terms of moments of the jump 
size intensity v. However, the jump intensity measure and the jump size intensity measure of 
the impulsive cylindrical process are closely connected, see Remark [TJ In particular, the jump 
intensity measure of the driving process (Zt)te[o,T] 111 ® i s allowed to be infinite. Thirdly, 
having found a suitable error expansion, the following estimates contain additional difficulties 
due to the infinite-dimensionality of our problem. This is why in Section 6 the assumptions 
on the jump intensity measure of the driving process are stronger than those in |28j . 
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2 Impulsive cylindrical process 

Let O c R d be open and bounded, T 6 (0, oo), and consider the product space [0,T] x O x R 
equipped with the Borel c-algebra i3([0,T] x x R). The generic element in [0, oo) xOxfi 
is denoted by (t,£,a) or (s,£, <r). Let v be a sigma-finite measure on R and 7r a Poisson 
random measure on [0, oo) x O x R with reference measure dtd£ v(da). By tt we denote the 
compensated Poisson random measure, i.e. 

■ft(dt, d£, da) = ir(dt, d£, da) - dtd^ v(da). 

Let 7r be defined on a complete probability space (f2,.A,P) with a filtration (J~t)te[o,T] satis- 
fying the usual hypotheses, cf. |23j . We assume that 7r([0,t] x B) is ( Tt )-measurable for all 
t e [0,T], B e i3(C x R) and that vr((s,t] x B) is independent of .F s for < s < t < T, Be 
B((D x R). Furthermore, we assume that tt is of the form 

oo 

7r(A)(w) = 2«(T i ( w ),B j ( w ),E J -( w ))(^) ) i^([O,T]x0xR) )We fl, (6) 

3=1 

for a properly chosen sequence {(Tj, Sj, Ej)) ^ of random elements in [0, T] x x R, cf. |27j . 
Chapter 6. 

For general impulsive cylindrical processes, v is a Levy measure, i.e. ^({0}) = and 
min (c^, l) v{da ) < oo, cf. [27J, Ch. 7.2. To simplify the exposition we will additionally 
assume that yi a .| >1 a 2 u(da) < oo which is equivalent to saying that 



/ a 2 v(da) < oo. (7) 



In the Appendix [Bj we show how our results can be extended to the general case. Under 
condition ([7J the random variables and Zt defined below have finite second moments. 

(This is, e.g. the case if v is a Levy measure with bounded support.) 

To fix notation, let us give a brief survey of L 2 -integration w.r.t. tt for deterministic inte- 
grands. Let / : [0, T] x O x R R be a simple function, i.e. 

n 

fc=l 

where n e M, a fc 6 R and A fe e <B([0,T] x O x R), such that 

[ T f f t Ak (t,£,a)v(da)d£dt<oo, 
Jo Jo Ju 



for k = 1, . . . ,n. 

For simple functions / the stochastic integral w.r.t. 7r is defined by 



Since 



Jr i r 11 
/ /(*, £, o") jr(dt, df , da) := V a fc 7r(A fc ). 

E lf T f f(t,t,a)TT(dt,dt,da)\ = [ T [ \f(t,t,a)\ 2 dtd£v(da), 

\J0 JOxR | JO JOxR 
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the stochastic integral can be uniquely extended to an isometric linear operator mapping 
L 2 ([0,T] x O x R, dsd^u(da); R) to L 2 (f2,.A,P;R). Later, we will also consider 7r-integrals of 
.ff-valued square integrable functions which are defined similarly. In the course of the proof 
of our result, we will have to deal with /^-integrals of stochastic integrands against the (not 
compensated) random measure 7r; see |27|, Sections 6.2, 8.7 and [33J, Chapter 2, for a detailed 
exposition of the various integrals. 

As usual, for t € [0,T] and / e L 2 ([0,T] x O x R, dsd^u(dcr);'R), we define 

f f f(s,£,a)n(ds,d£,da):= f f lr 0t i(s)f(s,£,a)7t(ds,d£,da). 
JO JOxE Jo JOxR L 1 

Now we can to define the impulsive cylindrical process. Let (e k )k£N be an orthonormal 
basis of H and for feel, te [0,T], let 

Z { t k) := Z(t,e k ) := /" f e k ((,) ■ afc(ds,d^,da). 
JO JOxR 

Note that the processes )te[o,T]> ^ e IN, are real- valued, square integrable Levy processes 
which are also martingales. 

Now let U be a further Hilbert space such that if is densely embedded in U and such that 
the embedding is Hilbert-Schmidt, e.g. U = H~2~ e ((D) for some e > 0. In |27| . Chapter 7, it is 
shown that 

n 

Z t :=L 2 (n,A,-P;U)-\im ^ Z\ k) e k , t€ [0,T], 
nt~ fc=1 

defines an [/-valued L 2 (P)-Levy-martingale with reproducing kernel Hilbert space ( ■ , ■ }%) = 
(H,J u a 2 u(da)(-,-) H ), H = L 2 (0). (Z t ) te[a T ] is called impulsive cylindrical process on 
L 2 (0) with jump size intensity v. 

Remark 1. The jump intensity measure \i of the [/-valued Levy-process (Zt) t€ r 0T i is given 
by 

fi(B) = \®v(u£,a)eOxR-. £e k (t)ve k 6 b\Y B € B(f7) 

where A denotes Lebesgue measure on O and the infinite sum is a limit in L 2 (OxR, d^u(da); U). 
As for the following examples, compare [25], [26J. 



Example 2. Let j3 e (0,2), r e (0, oo) and consider the jump size intensity 

1 
a~ 

Then, for every B e B(O), the one-dimensional process 

-t 



(z(t,t B )) u[0T] 



ts[Q,T] 



can be characterized in terms of one-sided /3-stable Levy processes as follows: 

If P 6 (1, 2), then (^(t, ls)) t£ r T j is equivalent to the process obtained by removing all jumps 

greater than r from a one-sided /3-stable Levy process (£t)t € [o,T] with Laplace transform 



Ee 



- rLt = exp {-t\B\c0 ■ r 13 } = exp Lt\B\ J™ (l - e'™ - ra)-^ J , r > 0, 



5 



(8) 



and then adding the shift t i-> t\B\f™ cr^f^ = Here \B\ denotes the Lebesgue measure 

of B e B(O). 

If (3 e (0,1), (Z(t, H-b))^^ T j is equivalent to the process obtained by removing all jumps 

greater than r from a positive one-sided /3-stable Levy process (£t)te[o,T] with Laplace trans- 
form 

Ee- rLt = exp {-t\B\c p ■ r?} = exp {-t|fl| £°° (l - e~ ra )-^ j , r > 0, 

and then subtracting the shift t >-> E Jq Jq J^^b(0 ' CT ^(ds,d^,d(j) = t\B\^ 
Example 3. Let a e (0, 2), r e (0, oo) and set 

^t* 7 ) = uiS 1 [-r,r]( ff )^ 

Then, for i? e i3(C), we have 

Z(t,l B )= [* [ t B (0-vHds,dZ,da) 

= L 2 (n,A,F; R)-lim /" f l B (f) ■ a ir(ds,d(,,d(r), 

e\0J0 JOx{\a\>t} 

and lB)) te j- T j is equivalent to the process obtained by removing all jumps of absolute 

value greater than r from a symmetric a-stable Levy process (-^t)t £ [o,T] with Fourier transform 

Ee irLt =exp{-t\B\C a -r a } = exp\-t\B\ /" (1 - cos(ro-))-^- 1 , reR, 

[ jk l°1 J 

i.e. 

Ee ^(t,l B ) =exp |_ t |_ B | ^ T (i_ cos ( r(J ))^_J ) r€E . 
The second equality in ([8| holds since f is symmetric and therefore 

f f 1 B (0 -adsdiv{da) = 0, e > 0, 

JO -/Ox{|o-|>e} 

so that the integrals w.r.t. tt and 7r coincide for integrands in L 2 ([O,T]x0xR, dsd£v(da); R)n 
^([OjT] x O x R, dsd^u(da); R). 

For the estimate of the weak order of approximation (Theorem |8|), we have to make a 
further restriction on the jump size intensity v by assuming that 

/ max (|o"|, <t 2 ) v(do) < oo. (9) 

This is, e.g. fulfilled if f3 e (0, 1) and a e (0, 1) in Examples 1 and 2. 
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3 Notation, assumptions and preliminary results 



The inner product and norm of H - L 2 (0) are denoted by (-,-}h an d | 1 Iff, respectively; 
Q : H -*■ H is a bounded nonnegative definite symmetric operator and A : D(A) c H H a 
(unbounded) strictly positive definite self-adjoint operator whose domain D(A) (endowed with 
the graph norm | • \h + \A ■ \h) is compactly embedded in H. Therefore, the spectrum consists 
of a sequence (\k)ktN c (0, oo)- We assume that the eigenvalues are ordered increasingly, 
Ai < A2 < • • •, including multiplicities. By (e^)/^ we denote the corresponding orthonormal 
basis of eigenvectors. For any s > we set 

{00 00 
u=Y.( u ^k)He k tH : Y, X k S (u,e k ) 2 H <ool, 
k=l k=l J 

00 

A s u := Y. X k{ u ^k)He k , ueD(A s ), 

k=l 

so that D(A S ) endowed with the graph norm | ■ \h + \A S ■ \h is a Hilbert space. Furthermore 
we define D{A~ S ) for s > as the completion of H with respect to the norm | • \h + \A~ S ■ \h, 
defined on H by \A- s u\ 2 H = ^^=1 K 2s (u,e k ) 2 H - 

Let (Vh)h>o be a family of finite dimensional subspaces of V ■= D(A 1 I 2 ) parameterized 
by a small parameter h > 0, given by the finite element method. As standard references for 
the finite element method we mention [5] and |31| . By Ph we denote the orthogonal projector 
from H onto Vh with respect to the inner product ( • , • )h and is the orthogonal projector 
from H onto Vh with respect to the inner product {A 1 ! 2 ■ ,A l l 2 ■ )h- For any h > 0, we define a 
positive symmetric bounded linear operator Ah ■ Vh -*■ Vh by 

(A h Uh,v h )H = (A 1/2 u h ,A 1/2 v h )H V(u h ,v h ) e V h xV h . 

As in [3, let (S(t))t>o denote the Co-semigroup of contractive operators on H generated by 
-A and let (Sh(t))t>o denote the semigroup of operators on Vh generated by -Ah- 

We assume that the finite dimensional spaces Vh admit the following two estimates, which 
are classical properties for standard finite element spaces. 

Assumption: 

For all q e [0, 2] there exist constants k\ > 0, k% > independent from h such that for all t > 

\\S h (t)P h -S(t)\\ L(H) < Kl hH^ 2 , (10) 
\\S h (t)P h -S(t)\\ L{HMAl/2)) < K 2 ht-\ (11) 

In [32], Theorem 3.5, [3], Theorem 3.2, and |18) . Theorem 4.1, the above estimates (|10p and 
(TTTj) are shown for the case (A,D(A)) - (-A, H 2 {0) n Hq(O)) under the assumption that 

\U h v-v\ H < K h s \A s / 2 v\ H , (12) 
\A l ' 2 (Ti h v-v)\ H < Koh^A^vlH (13) 

for all s € [1,2], v e D(A S I 2 ) and some constant kq > 0. Finite elements satisfying (jl 2|) and 
(|13p are well known, like P k triangular finite elements on a convex polygonal domain or Q k 
rectangular finite elements on a rectangular domain, k > 1, see |5], [31] . 
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The main assumptions concerning the operator A and the covariance operator Q of the 
noise are the following: There exist real numbers 

a>0, /3e(a-l,a] (14) 

such that 

oo 

Tr(A-°) = Z K a < oo , (15) 

71=1 

APQtL(H). (16) 

Notice that (|16p implies for any A e [0, 1] 

A X/3 Q X e L(H) and \\A X ^Q X \\ L(H) < \\A^Q\\ X L(H) . (17) 

According to [27j, Chapter 9, equation Q has a unique (predictable) weak solution which is 
given as the mild solution 

X t = S(t)x + £ S(t-s)Q 1/2 dZ s , te[0,T], (18) 

provided that 

l|5(t)Q 1/2 ||(HS)eL 2 ([0,T],dt), (19) 

where || ■ || (hs) ^ s the Hilbert-Schmidt norm. It is shown in [9] that f|14J> . (j!5p and (|16jl are 
sufficient conditions for (|19p . Throughout this paper, we only consider weak respectively mild 
solutions of SPDEs. 

We proceed with some further notation: By Cj;(H) = C 6 fc (if;R) we denote the space of 
all fc-times continuously Frechet-differentiable real valued functions on H which are bounded 
together with their derivatives. For (j) e Ch(H) and x e H the first-order derivative D<p(x) 
of (ft in x is identified with its gradient and thus considered as an element in H. Similarly, 
for <ft e C 2 (H) and x 6 H, the second-order derivative D 2 tp(x) is seen as an element in 
L(H) = L(H,H), the space of all linear and bounded operators on H. 

Let L( HS )(ff) = L(ft S }(H, H) denote the Hilbert space of all Hilbert-Schmidt operators on H. 
It is a subspace of L(H). Given an orthonormal basis (ek)kaN of H, the scalar product in 
L (HS)( H ) °f operators T e L (HS) (#), S € L (HS) (F) is given by 

( t ; 5')(hs) : = E^e^Te^. 

fceW 



Finally, we use .My(-ff) for the space of all right continuous L 2 (P)-martingales M = 



The corresponding Hilbert-Schmidt norm is denoted by || ■ || (HS) 

(Mt)te[o,T] with values in H. ([^f]t)te[o,T] denotes the tensor quadratic variation (or operator 
square bracket) of Me .M^-ff); ([-^]t)t 6 [o,T] denotes the continuous part. Note that for P- 
almost every u e 0, for every £ 6 [0,T], [Mjt(u;) and [Mj£(a;) belong to the space of nuclear 
operators on H, a subspace of L/gg)(iT) which is continuously embedded in L^^^H), cf. |23j . 

4 Approximation scheme 

In order to approximate the mild solution 

X t = S(t)x + f * S(t-s)Q 1/2 dZ s , t€[0,T], 
Jo 
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of equation ((4]), we adapt the numerical scheme described in [9], with (V^t)t 6 [o,T] replaced by 
the jump process {Z t ) t e[o.T]- 

Given an integer N > 1, set At = T/N and t n = nAt, n = 0,...,N. For any h > 0, the 
approximations X£ of X tn in Vh, n = 0, . . . , N, are defined as those V/j-valued random variables 
which satisfy for all Vh e Vh 

{xr 1 - XI , v h ) H + At [A^{6Xt l + (1 - 0)Jtf ) , ^V)^ = (Q 1/2 Z tn+1 - Q l ' 2 Z tn , v h ) H 

(20) 

and 

(Xlv h ) H = (x ,v h ) H , (21) 

with 

9e (1/2,1]. (22) 
Here the expression (Q 1 ^ 2 Zt n+1 - Q 1 l 2 Z tri , Vh) H is the usual shorthand for 

M 

L 2 (n,A,P;R)- lim £ (Z(^ + i,e fe ) - Z(t n ,e fe )) • Q 1/2 e fc , v fc L , (23) 

( e fc)fc£M being an arbitrary orthonormal basis of H. In particular, one can choose (ek)k<=ff as 
an orthonormal basis consisting solely of eigenvectors of Q and of elements of the kernel of Q. 
Note that (XJ^) ne ^ ^ N y is a discretization of (^t)t<=[o,T]> both in time and space. Solving the 
equations above leads to 

n-l 

K = St At P h x Q +Y.SlM X ThAtPhQ X, \Z tk+x -Z th ), (24) 

k=0 

for 7i € {0, ... , iV} and /i > 0, where 

Sfc,A t := (/ + ^AtA)^(/-(l-^)AtA), 

r h ,At == (/ + 0Aa h ) _1 . 

Remark 4. In this paper, we are interested in the weak order of the scheme ()20p -(|2i p . Theo- 
rem [8]below gives an estimate for the difference \Eiip(X^) - Wj(p(Xx)\ for suitable real valued 
functions tp. Obviously, for the numerical implementation one has to truncate the infinite sum 
in (|23p . That is, one has to replace the right hand side of (|20p by 

({Q^) 1/2 Z tn+1 - {QM) 1/2 Z tn , v h ) = £ {Z(t n+1 ,e k ) - Z(t n , e k )) ■ (Q^e k , « fc ) , 

k=l 

for some Mel, Here, the operator Q^ M ^ is defined by 

M oo 

fe=l fc=i 
Let fxr , '^ M ^) » rl be the solution of the truncated scheme, i.e. 

V fi /ne{0,...,A r } ' 

n-l . 

K (M) = si At p hX0 + ^ ^r^Y^i - z ^ 
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and let 

X { t M) = S(t)x + £ S(t - s)(Q^) 1/2 dZ s , t e [0, T]. 
Then, in the setting of Theorem [8] below, it is not hard to see that 

|E^(xf' (M) ) - Ey?(4 M) )| < C ■ (hP + (At)^), 

where C > and 7 > are the same numbers that appear in the upper bound for the error 
\E<p(Xfr) -~E<p(Xt)\ in Theorem [8] In particular, the constant C does not depend on M 6 IN. 

Therefore, the error \~E<p(xj*'^ M ^ — Ey^-Xr)), which is relevant for numerical simulations, can 
be estimated by 



\~E<p(x" ,CM) ) - J£<p(X T ) < C • (/1 27 + (At) 7 ) + Ev?(X^ M) ) - E^(X T ) 

It will be useful to introduce some further notation. For each h > 0, we will also consider 
the following spatial discretization of (Xt)tz[o,T] 

X htt = S h (t)P h x 0+ [ l S h {t-s)P h Q 1 l 2 dZ s , ie[0,T], (25) 
Jo 

and the two auxiliary processes 

Y h)t = S h (T)P h x + £ S h {T - s)P h Q l l 2 dZ s 

= S h (T)P h x + f ${s)dZ s , te[0,T], (26) 

rt N-l 

Yh,t = S% At P h x + £ Slti 1 T hA tMt k ,t k+1 ]{s)P h Q ll2 dZ s 
J0 fc=0 

= S% At P h x + J*T(s)dZ s , te[0,T]. (27) 



To ease notation, we have used 



:= S h {T-s)P h Q 1 l 2 , 

N-l 

r(*) == Y, S hAt~ lT hAt 1 (t k ,t k+1 ] (s)P h Q lf2 . 

k=0 

Moreover, it will be convenient to set 

$(s) :=*(T-a). 

Before we come to the results, we let us rewrite the stochastic integrals in a way that fits 
to our purposes. The integrals with respect to dZ s can be written as 7r-integrals. To this end, 
let us define mappings E, F, G, F from [0,T] x0xfi into H, by 

00 

E(s,Z,a) == Y, e k(0<rS(T-s)Q 1 ' 2 e k , (28) 
k=l 
00 

^0,^) == Z e k(0<rHs)e k , (29) 
fe=l 
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G(s,£,a) := Y, e k(Oo-T(s)e k , (30) 



k=i 



F(s^,a) := Y, e k(0^(s)e k , (31) 

k=l 

where (ek) k ^ is an orthonormal basis of H and the infinite sums are limits in the space 
L 2 ([0,T] x O x R, dsd£v(da); H). These limits exist since the operator valued functions 
s h> S(s)Q 1 / 2 , s *-> <3?(s) and s T(s) belong to L 2 ([0,T], ds; L( HS )(ff)) and because of the 
integrability assumption (jTJ). As usual, we do not distinguish between measurable mappings 
[O,T]x0xR H and their ds d^ z;((ier)-equivalence classes. Notice that the following equalities 
hold: 

f t S(T - s)Q l/2 dZ s = f f E(s,Z,a)Tt(ds,dt,da), 

J0 J0 JOxR 

f $(s)dZ s = f f F(s^,a)7t(ds,d(,da), 

JO Jo JOxR 

[ t T{s)dZ s = r f G(a,Z,a)*(d8,dt,da), 
Jo Jo JOxR 

f $>(s)dZ s = f f F(a,Z,a)Mda t d£,da). 

JO Jo JOxR 

Indeed, if we consider for example Jq &(s) dZ s , it is easy to verify that 

[ t ^{s)dZ s = L 2 (Q,A,F;H)- lim V f* $(s)e k dZ^ k) 
Jo n^'x'jZ^Jo 

n rt r 

= L 2 (n,A,¥;H)- lim V / / e fc (f )<T$(s)e fc 7r(ds,d£,d<r) 

= f f I E e fc(0^^( s )efe| Tr(ds,dtda), 
Jo JOxu\^[ J 

where the infinite sum in the last integral is a limit in L 2 ([0,T] x O x R, dsd£v(da); H). 

5 Error expansion 

In this section, for suitable functions (p defined on H, we state and prove representation 
formulas for the time discretization error Ey?(A^) - Eip(X h:T ) and the space discretization 
error Ec^A^^-Ec^A^). The errors are represented in terms of the functions v% ■ [0, T]xH 
R and v h : [0, T] x H -+ R defined by 

:= Ec^x + J^S^T-rOP^Q 1 / 2 ^), (32) 

u(i,x) := E<p[x + jJ t S(T-r)Q 1/2 dZ r y (33) 

For the reader's convenience, we summarize all equations, definitions and assumptions which 
will be needed in the formulation of our error representation theorem below: 
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Summary / Assumptions: 

• The equation we are interested in is 

dX t + AX t dt = Q 1/2 dZ t , X = x £H, te[0,T]. © 

• Assumption: The jump size intensity v of (-^t)t £ [o,T] satisfies 

f a 2 u{da) < oo. © 

• Assumption: The mild solution of equation (|4j) exists, i.e. 

iis(t)Q i/2 ii ( HS)^ 2 ([o,n^) <m 

(This ensures that the weak solution exists and is equal to the mild solution.) 

• We consider the following space-time- and space-discretizations of the weak solution of 
® 

n-l 

n = S n hAt P h x + Y, Sl% l T hAt P h Q l l 2 {Z tk+1 - Z th ), n e {0, ... ,7V}, (JMD 
X h!t = S h (t)P h x + £ S h (t-s)P h Q 1 ' 2 dZ s , ie[0,T]. (J2SJ) 

Theorem 5. Assume that and (fT9|) /io/d. Let ip e C 2 (H) such that for each x e PL the 
derivative D 2 (p(x) is an element of L( H gj(H) and that the mapping x h> D 2 tp(x) e L{- HS j(H) 
is uniformly continuous on any bounded subset of PL . Let T > 1 and (Xt)te[o,T] be the PL-valued 
stochastic process satisfying equation (j4|). For any N > 1 and h > 0, let (X^) ne ^ N y be given 
by < \24\ ) and let {Xh,t)u[o,T] be as in i\25\) . 
Then the following error expansions hold: 

Kip(X^) - E V {X h)T ) = {v h (T, S% At P h x ) - v h (T, S h (T)P h x )} 

+ vfj {v h (T - t, Y h>t - + G(t, £,a))- v h (T - t, Y h>P . + F(t, £, a)) 

+ (D x v h (T - t, Y hjt _), F(t, £, a) - G(t, £, aj) R } dtd£ u{da) 
=:/ + //, (34) 

^{X KT ) - E<p(X T ) = {v(T, S h (T)P h x ) - v(T, S(T)x )} 

+ ^fj {v(T - t, Y h>t - + F(t, C, a)) - v(T - 1, Y M _ + E(t, £, a)) 

+ (D x v(T - t, F M _), E(t, e, a) - F(t, £, <r)) ^ } dt d£ u(da) 

=: /// + LV. (35) 

Remark 6. In the following proof of Theorem [5j the uniform continuity of the mapping 
H 9 x h> D 2 (p(x) e L(HS)(-ff) on bounded subsets of H is needed to be able to apply Ito's 
formula in infinite dimensions, see [23], Section 27. This assumption is always fulfilled in finite 
dimensions and in literature it is sometimes forgotten to mention in the infinite-dimensional 
case. 
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Proof of Theorem^ We begin with the time discretization error Eip(X^) - E^X^.t)- Due 
to the definition of we have 

Mp(X h , T ) = v h (T,S h (T)P h x ), (36) 
Eip(X^) = Eip(Y hiT ) = Ev h (0,Y hjT ). (37) 

Applying Ito's formula to the function (t,x) Vh(T - t,x) and the (cadlag) martingale 
( Y h,t)te[o,T] yi eld s 

v h (0,Y htT ) = v h (T,Y h!O )- / —j-(T-t,Y ht _)dt 

Jo ot 

+ fj {D x v h (T - 1, Y h>t _), dY h ,t) H +UJ {D 2 x v h (T - 1, Y hjt _), d{Y h) .^) m 
+ £ K( T - t, Y h , t ) - v h (T - 1, F M _) - (D x vh(T - i, F M _), AY M ) ff }. 

(38) 

A process is called cadlag (continu a droite et limites a gauche) if almost every path is right- 
continuous and has finite left limits. We have used the standard notation Y/jj_ = lim s /tY^^ 
and AYh t t = Y^t - Yhj- which make sense for every cadlag process. We will now consider the 
mean values of the terms on the right hand side of (|38p separately: 

(i) Ev h {T,Y hjQ ) = Ev h (T,S% At P h x Q ) = v h (T,S% At P h x ). 

(ii) E / T (D x v h (T - t,Y h}t _), dY Kt ) H = because (f* (D x v h (T - s,Y hjS _), dY h,s) H ) te[QT] is 
a martingale starting in 0. 

(iii) ([P\.lt) t6 rQ T j = because (?/i,t)fc[o,T] i s a purely discontinuous martingale. This follows 

because {Yh,t)u{o,T] ls tli e Ai\{H)A.\mit of finite sums of its coordinate processes which 
can be easily identified as purely discontinuous martingales. Since the space of all purely 
discontinuous L 2 (P)-martingales A4^ d (H) is closed in J\4^(H), the process (Yi,t)te[o,T] 
is also purely discontinuous (cf. |23| . Chapter 4). Hence ([[Y/^.Jj ) U [ T ^ = 0j an d therefore 

±Ef o T (D 2 x v h (T-t,Y h ^), 4n,]?) (HS) =0. 

(iv) Concerning the jump term 

E Y, {v h (T - t, Y h , t ) - v h (T - 1, Y M _) - {D x v h (T - t, Y M _), AY h)t ) } (39) 

t<T 

we give only a brief sketch how to deal with it; the full argument is postponed to the 
appendix. 

Let m be the jump counting measure of (Yh t t)te[o,T]i i- e - 

s<t 

The basic idea is to write the sum in (|39p as a (pathwise) integral with respect to m and 
then to consider the random measure m as an image measure of tt with respect to the 
function G defined by (|30|) : 

E £ {v h (T - t, Y Kt ) - v h (T - 1, Y M _) - {D x v h (T - t, Y M _), &Y U ) R } 

t<T 
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-T 

= E 



-T 

E 



f jf {«h(T - t, + y) - ifc(T - t, F M _) - (D x v h (T - t, Y M _), } m(dt, dy) 

T T £ {« fc (r-t,y h)t _ + G(t,e,<r))-t;fc(r-t,y Ait _) 

JO JOxR «■ v 7 

-(A,i>h(r-t,Y hjt _), G(t,e,a)) H }vr( a !t,de,dc7). (40) 

Here the last term is the expectation of an stochastic L -integral with respect to the 
random measure ir, cf. |27j . Section 8.7. Since dtd^i>(d<j) is the compensator of 7r, 

e( T ( W(r-t,y M _ + G(i,e,a))-^(T-i,r M _) 

JO JOxR I V 7 

- {D x v h (T - t, Yh,t-)i G(t,t, a)) H } ir(dt, d£, da) 
- IE f T f {v h (T-t,Y h)t _ + G(t,Z,<r))-v h (T-t,Y hit _) 

JO JOxR l V 7 

-(D x v h (T-t,Y h , t -), G(t, a)) H )dt divider). 

(v) Going back to the usual construction of the stochastic integral, one can see that the laws 
of the random variables 

fJ t S h (T-s)P h Q 1 / 2 dZ s and f* S h (s)P h Q 1/2 dZ s 

are equal. Consequently 

v h (t,x) = Ecp(x + £ S h {s)P h Q l l 2 dZ s \ 

= E<p(x + J t ^{s)dZ^, 

where we have used the notation of Section H $(s) = $(T-s), $(s) = S h (T - s)P h Q 1/2 
for s e [0,T]. Now we apply the ltd formula to the function H 9 y h> + y) e IR and 
the (cadlag) martingale $(r) dZj.)^^ ^ 

4rKhx r ^ (s) ^)' e(r) "«> (HS ) 

+ £ ( x + X" dZs ) ' ip { x + £ rfZs ) 

-[pv(x + J* Hs)dZ,y a£ ^dZ^J. 

Here, / Q r $(s)dZ<j denotes the (pathwise) left limit lim q/ , r $>(s) dZ s and the expres- 
sion ^(r)® 2 in the second line of the formula stands for the operator L^^(H) -*■ 
L(hs)(-^)) T h- > <I>(r)T<I>(r)*. Taking the expectation on both sides of the equation, 
reasoning as in (iv) and differentiating with respect to t yields 

d r 

—v h (t,x) = / {v h (t,x + F(t,£,a))-Vh(t,x) 
at JOxR 

- {D x v h (t, x), F(t, e, a)) H } di u{da). (41) 
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Finally we get flMJ) if we combine (J36J) - ([38]) and (i) - (v). 

Now consider the spatial error E^X^^-E^X^). The definition of the function v implies 

E<p(X T ) = v(T,S(T)x ), 
E<p(X h , T ) = Ev(0,Y KT ). 
As above, we can apply Ito's formula to (t,x) *-> v(T-t,x) and (Yh t t)te[0,T] to get 

~ T dv 



r 1 ov 

v(0,Y htT ) = v(T,Y h>0 )- / — (T-t,Y hyt _)dt 

Jo at 

+ J o (D x v(T-t,Y h>t _),dY h>t ) H 

+ \£{Dlv{T-t,Y h ^),d{Y h X) {nS) 

+ Y, \v(T - 1, Y h>t ) - v{T - t, y M _) - (D x v(T - 1, AY h>t ) H }. 

t<T 

This can be used, similar to the argumentation leading to (|34p . to verify (1351) . □ 



Remark 7. The proof of Theorem [5] reveals that we can extend Theorem [5] to equations of 
the form 

dX t + AX t dt = Ql /2 dW t + Q{ /2 dZ t , X -xq€ H, te[0,T], (42) 

where (Wt)ts[o,T] * s a cylindrical Wiener process on H which is independent of (^t)t e [o,T] 
and the covariance operators Qq and Q\ are, just as Q above, bounded, nonnegative definite, 
symmetric and satisfy (|14p . (|16p . The solution (Xt)t s [o,T] °f ( BZD is given by 

X t = S(t)xo + £* S(t - s)Q{ /2 dW a + £* S(t - s)Ql /2 dZ s , te[0,T], 
and the discretizations are 

n-l 

K = Sl M P h x +Y s h;M l ThAtPhQT(W tk+1 -W tk ) 

n-l 

+ £ ShAt^hAtPhQiiZt^ - Z tk ), h>0,ne{0,...,N}, 

k=0 

X h , t = S h (t)P h x + £* S h (t-s)P h Ql /2 dW s 

+ £ S h (t-s)P h Q{ /2 dZ s , te[0,T]. 

Assuming for simplicity that Q\ = Q2, the time discretization error ~Eip(X^) -~E(p(Xh,T) anci 
the spatial error 'Eip(XhT) ~ E^Xy) are obtained by adding the terms 

^E Jj Tr{r(i)* D 2 x v h (T -t,Y h ^)T(t) - *(*)* D 2 x v h (T - t,F M _)$(t)} dt (43) 

and 

~E Jj Tr{<&(i)* D 2 x v(T-t,Y h ^)$(t) - (S(T -t)Q l ' 2 y ' D 2 x v{T -t,Y h>t _)(S(T -t)Q{ /2 )}dt 

(44) 

to the right hand side of ([34 1) and (135 1) . respectively. (Here one has to replace Q by Q\ in the 
definition of $(t) and T{t).) 
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6 Weak order of convergence 

In this section, we show an estimate for the weak order of the convergence of to X? as h 
tends to zero and N tends to infinity, given the integrability condition (|9]) . The proof is based 
on the error expansions in the last section. 

For the convenience of the readers we summarize the assumptions made in the theorem below. 

Assumptions: 

• The jump size intensity v of (-^)te[o,T] satisfies 

max (|cr|, a 2 ) v(do) < oo. ([9]) 

• The finite element spaces Vh, h > 0, are such that for all q e [0,2] there exist constants 
K\ > 0, k-2 > independent from h such that for all t > 

\\s h (t)p h - s(t)\\ L(H) < K X hH-*' 2 , m 

\\s h {t)p h -s{t)\\ L {.HMA^))^^ hrl - (ED 

• There exist real numbers a and /3 with 

a>0, /3e(a-l,a], (JUJ) 

oo 

Tr(A~ a ) = £ K a < oo , (USD 

71=1 

aPq£L(h). (HH) 

(Note that this particularly implies assumption (jl 9[) needed for Theorem [5]) 

Theorem 8. Assume ©, (fT0]t . (jTTj) . (fT4|) . (fT5|) and (Jl6j) feied a&oue. Lei 99 e C 2 (F) suc/i 
t/ioi for each x € H the derivative D 2 cp(x) is an element of L( H gj(H) and that the mapping 
x h> D 2 ip(x) e L(Hgj(H) is uniformly continuous on any bounded subset of H. Let T > 1 and 
(^)ts[0,T] ^ e ^ e H-valued stochastic process satisfying the equation 

dX t + AX t dt = Q 1/2 dZ t , X = x €H, te [0, T] . © 

For any N > 1 and h > 0, Zef (^^) n e{o,...,Af} ^ e 9"i ven by (\24\)- 

TTien /or any 7 < 1 - a + < 1, there exists a constant C = C(T,ip,A,Q,\xq\,'),v) > which 
does not depend on h and N such that the following inequality holds 

\mp{xj?) - Ev9(x T )| < c ■ (h 2 -< + (Aty), 

where At = T/N < 1. 

Proof. As in Theorem [5j we split the error into the time discretization error and the spatial 
error, 

mp(X% ) - Ecp(X T ) = {E^(Xf ) - E<p(X h , T )} + {Ev?(AV) - Eip(X T )}. 

We will estimate each term separately. Throughout this proof, C denotes a positive (and finite) 
constant that may change from line to line. 
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Consider the time discretization error given by ( I34|) . 'Eip(Xj^) -~Eip(Xh^) = I + 11. Clearly, 
due to the definition of Vh, 



\ J \ * h\\cl(H) \\ S h,At P h ~ S h(T)P h \ L{H) \x Q \ H . 

Using spectral calculus and (|22l) . \Sf^^ t Ph - Sh(T)P}^ can be bounded uniformly with 
respect to h. To be more precise, for T > 1 we have 



S? At P h -S h (T)P h \ 



T(m < sup 

L{H > A>0 



sup 

r>0 



-NXAt 



/ l-(l-g)AAt \ 

\ i + e\At ) 

N 



N 



_ Nr / l-(l-fl)r \ 
\ l + Or J 



C 



< — < C ■ At, 

N 

cf. [20] . p. 921, Theorem 1.1 for the penultimate estimate. One obtains 

|/| <C- At. 



(45) 



Concerning II, we apply the mean value theorem and the Cauchy-Schwarz inequality and 
obtain 

II < f T f 2\\D x v h \\ Cb(H) \F(t,£,tT) - G(t,t,a)\ dtd^u(da) 



nD x v h \\ Cb(H) f R \a\,(da) fj f o 



Ee*(0(*(*)-r(t))e* 



k=l 



dtd£, (46) 



H 



where the expression E^Li efc(£) (<&(<;) -T(t))efc has to be understood as the value at (t,£) of 
a (it cZ£- version of 



/ N 

L 2 ([0,T]xO,dtdH; H)-lim (t,£) " £ e *(0($(*) " 

JV - 00 V fe=i 

Next, let 7 > and 71 > such that 0<7<7i<l-a + /3<l. Using 
00 

£e fc (0(<&(t)-r(t))e fc = 



fc=i 



JV-l 



n=0 

6|) can be estimated by 



^(T-i)-£l (Wn+l] (i)S£^^ 

k=l 



N-l 



II <C I [ Sh (T-t)-Z t {tnU+l] {t)S^ t - l T hAt \A' h 



n=0 



.(l-7i)/2 



dt 



L(V h ) 



O 



k=l 



d£ 



H 



<c 



I, 



T 



N-l 



s h (T-t)~ £ l^wiWflKr 1 ^* 4 

n=0 



(1-71 )/2 



(ft 
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E e t (o4 71 " 1)/2 fl.e l/2 « 

fe=i 



1/2 



= C\\A^I 2 P h Q^\\ 



(HS) 



r 



N-l 



S h (T-t)- £ l^t^S^T^AA 



n=0 



|(l-7l)/2 



(ft. (47) 



Following E|, ||4 71 " 1)/2 PhQ 1/2 L HS) 

can be bounded from above by some constant C > which 
does not depend on the choice of h > 0: 



\A ( r m r h Q 1/2 \\ m ±c. 



(48) 



(49) 



This is essentially due to (|15p and (j 1 6 1) . Furthermore, again following [9], 

(s,(r - o - i> (w „ +l] (i)^r%At) 4 1_7l)/2 

V n=0 / L(V h ) 

< | CAt^J^i 2 l (tlliWl] (t)((JV'-n-l)At)-(( 1 -^)^) teCO,^-!] 
" \ C(T-t)^/ 2 , t€(t N -!,T] 

which follows directly from spectral calculus. A combination of (|47p . (|48]) and (|49p yields 
^ [ E l (tn|t „ +l] (i)((^-^-l)At)- ((1 -' yi)/2+7) di+ / (T-i) (7l - 1)/2 ^ 

J0 n=0 ~>tN-i ) 

<C- At 7 , 

as (1 - 7 i)/2 + 7 e (0, 1) and (T - t)^- 1 )/ 2 < (T - t)-((i-7i)/2+7) f or t € (t^.^T], 
Finally JSJ), (^J and JSJJl add up to 



(50) 



|E^(X,f)-E^(AV)|<CAt 7 



(51) 



for all T > 1 and At < 1. 



Now we turn to the spatial error given by ([35 1) . Ey^X^r) - E(/j(Xr) = III + IV. Firstly, 
according to the definition of v, (|33[) . and resulting from (jlOp with g = 27 < 2, 



|///| < \\<p\\ cUH) \\S h (T)P h -S(T)\\ L(H) \x \ H 
^ WfWcUH) (Kih 2l T-~<) \x Q \ H = C-h 2 \ 



(52) 



Considering IV, we apply the mean value theorem and the Cauchy-Schwarz inequality and 
obtain 

" f 0xR 2\\D x v\\ Cb{H) \E(t,C,a)-F(t,C,a)\dtdCu(da) 
= 2\\D x v\\ Cb (H) I \o\v{da) 



/ Z e k(0{S(T-t)-S h (T-t)P h )Q l ' 2 e k 

JO k=l 



dtd£, 



(53) 



H 
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where the expression XifcLi e k(£){S(T - 1) - S^iT - t)P^Q 1 l 2 e} t has to be understood as the 
value at of a dt d£- version of 

L 2 ([0,T] x C,dtd(; H)- lim l(t,0 ~ £ e fc (£)(S(T- 1) - S h (T -t)P h )Q l ' 2 e k \ . 

N ^°°\ k =i I 

Pick 71 > such that 0<7<7i<l-a + /3<l. Because of ([17} we know that A^^Q 1 / 2 € L(H). 
Furthermore, the fact that 1 - 7i + /3 > a implies A~ (1 ~ 7l+/3)/2 e L {m) (H). It is also not hard 
to see that the operator (S(T - t) - Sh(T - t) P^) A^" 11 ^ 2 has a continuous extension defined 
on the whole space H. Therefore we may write 



£ e k (0(S(T -t)- S h (T - t)P h )Q^ 2 e k 

k=l 

oo 

= (5(T - 1) - S h {T - t)P h ) A^» 2 £ e fc (0^ (1 - 7l+/3)/2 ^ /2 Q 1/2 e fc . 

k=l 

A further application of the Cauchy-Schwarz inequality in (|53p gives 

- T 



1/2 



L{H) 



dt 



\IV\ < CJ Q \\(S(T-t)-S h (T-t)P h )A^/ 2 

x(Xmo^) 1/2 (/ c 

We are now going to show that for all t e [0, T] 
J (s(T-t) - S h (T-t)P h )A^' 2 \\ L(H <Ch 2 ~< (r( 7l(7l - 1)/(27)+1 ) + r( (1 - 7l)/2+7 )) . (55) 

This will be done in several steps. 

(i) Due to the self adjointness of S(T - t), S h (T - t)P h and A^ 1 ^ 2 , 



L(H) 

oo 

£e fc (0^ (1 - 7l+W2 ^ /2 Q 1/2 e fc 

11)12 
^> 2 dt. 

L(H) 



(54) 



\(S(T-t)-S h (T-t)P h )A^ 



-7l)/2 



L{H) 



{{S{T-t)-S h (T-t)P h )A^I 2 ) 
A^/ 2 (S(T -t) - S h (T -t)P h ) 



\L(H) 



L(H) 



Here the operator 

^(i-7i)/2(5( T _ t ) _ Sh (T _ t ) Ph ) is properly defined since S(t)H 
D(A X ), t > 0, A e R and V h c D(A 1 / 2 ) c D^ 1 " 71 ^ 2 ). 

(ii) Due to the Holder inequality we have for all x e D(A 1 I 2 ) 

\A^I 2 x\„<\A 1 l 2 xt 11 -1*12. 



H WH- 



Consequently, 

\\A^/ 2 (S(T -t) - S h (T -t)P h )\ 



L(H) 



< \\S(T - t) - S h (T - t)P h \ X -^ D{Am)) ■ \\S(T -t)- S h (T - t)P h \l {H) . (56) 
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(iii) If 27 + 71 > 1, then (27 + 71 - l)/7 e [0, 2]. In this case, we combine (|56p with (jlip and 
(TTU|) . q = (27 + 71 - l)/7- Provided that h < 1, one gets 



\L(H) 



\\A^-^ 2 (S(T-t)-S h (T-t)P h )\ 

< C ( h l ~~ n t 11 _1 ) ( /l^ 27+71 _1 )T 1 /T^( 2 T + T 1 -1 )ti/( 2 7) ^ 

< c7i 27 r( 7l(7l_1)/(27)+1 ). 

(iv) If 27 + 71 < 1, we first combine (|56p with (jlip and (|10p choosing g = and get 

|| A d-7i)/2 ( 5(r _ t) _ Sh{T _ t )P h )\\ L{H) < Ch 1 ^ (T - i)" (1 " 7l) . (57) 

Secondly, again using (|10p with q = 0, one derives 

|| A (l-7 1 )/2( 5(T _ t) _ 5/i(T _ W )|^ ) 

< \\S(T - 1) - S h (T - t)P h \ l L ^ D{A1/2)) \\S(T - 1) - S h (T - t)P h \\l\ H) 

< C\\S{T-t)-S h {T-t)P h \ l L ^ D{All2)) 

< C (\\A l l 2 S{T - t)\\ L{H) + \\A^ 2 S h (T - O^lk(H)) 1 " 71 
= C(\\A l l 2 S{T - + \\A]! 2 S h {T - O^IU^)) 1 " 71 

< C{T-ty^ l)l2 . (58) 

In the last step, the inequality 



supx e e tx <[-) t \ t > 0, e > 

x>0 



0' 



has been used. 

Now let A := 27/(1 - 71) 6 (0, 1). Interpolating J5ZJ) and (|5S]) yields 



1-A 



Combining (i), (ii) and (iii) yields (|55p . 

As 71(71 - l)/(27) + 1 < 1 and (1 - 7i)/2 + 7 < 1, we can integrate (|55p with respect to the 
time variable t. Then (15511 . (1521 . (l54|) and (1551) add up to 



|E^(X ft)T )-E^(X T )| < C/i 27 . 

□ 

Remark 9. (i) The proof of Theorem [8] can easily be combined with the proof in [9] in 
order to obtain the same result for the equation 

dX t + AX t dt = Q l J 2 dWt + Q^ 2 dZ t , X = x tH, te[0,T], (59) 
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where (Wt)t€[o,T] is a cylindrical Wiener process on H which is independent of {Zt)u[o,T] 
and where the covariance operators Qq and Qi are bounded, nonnegative definite, sym- 
metric and satisfy (|14p . (|16p . The corresponding discretization and error expansion has 
been described in Remark 1. 
Of course, the result also holds for the equation 

dX t + AXtdt = Q 1 /2 dW t + Q 1 1 /2 dZ 1<t + Q 1 2 /2 dZ 2>t , X = x €H, ie[0,T], (60) 

where (Zi,t)i 6 [o,T] an d (■^2,t)te[o,T] are impulsive cylindrical processes satisfying the 
condition required of (Z t ) t e[o,T] ( see ® above), the processes (Wt)te\oT]> (■^l,t)te[0,T]i 
(^2,t)t£[o,T] are independent, and the covariance operators Qq, Qi and Q2 are bounded, 
nonnegative definite, symmetric and satisfy (|14p . (|16p . For example, one could consider 
impulsive cylindrical processes described by the jump size intensities 

v%{da) = 1 l[_ T|r ](o-)d(T, i = l, 2, 

with indices of stability < ol\ < a.2 < 1. 

(ii) One might try to avoid the integrability assumption Q in the proof of Theorem [8] by 
rewriting the terms II and IV in (|34p and (|35p using Taylor's theorem. Obviously, 



II = E f f {(G(t,Z,a),D 2 x v h {T-t,Y h ^ + 9G(t,Z,a))G(t,Z,a)) 
JO JOxE 

- {F(t,£,a) , D 2 x v h (T-t,Y h;t - + OF(t,Z,a))F(t,S,a)) H }dtd£v(da), 

IV = E f T f {{F(t,Z,a),D 2 x v(T-t,Y h ^ + $F(t,Z,v))F(t,Z,a)) H 
JO JOxti 1 

^ t ,(r-i,Y h)t _ + ^(t,$,(r))S(t,$,ff)) H }dtdei/(d(r), 

where # = 0(t,^,u) e (0, 1) and # = i?(t,£,u;) e (0, 1). However, the integrands appearing 
here cannot be estimated analogously to the estimates of the integrands in the terms 
(|43[) and (|44p in which appear in the case of Gaussian noise. The reason is, that in 
the case of impulsive noise one has to estimate the difference of the second derivatives. 
As one needs integrability of suitable upper bounds for the integrands, this leads to 
profound technical complications. 



A Proof of equality ( 1401 ) 



Here we give a detailed proof for equality (|4U|) . We have to show 

E £ {MT - t, Y h , t ) - v h (T - 1, Y h ,t-) ~ {D x v h (T - t, Y M _), AY hjt ) } 

t<T 

= E f T f {v h (T-t,Y h:t _ + G(t,C,(r))-v h (T-t,Y htt _) 
JO JOxR I v ' ' 

- (D x v h (T - 1, G(t, M) H } n(dt, d£, da). 

We have already seen that 

Yh,t-Sh,te p hXQ= f T(s)dZ s = f f G(s,£,a)7r(ds,d£,dcr), 

JO JO jOxR 
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with G(s,^,a) - Efc°=i e k(0 <J ^( s ) e k^ where the the infinite sum is a limit in the space L 2 ([0, T]x 
O x R, ds d^u{da)\ H). By a standard monotone class argument, there exist simple functions 

m(n) 

G n (s,Z,a) = £ a n>k t An>h (s,Z,a), (s,£,a) e [0,T] x x R, 
fc=i 

where n e IN, m{n) e IN, a n) fc e and A nk e i3([0,T] xOxE) for k = 1, .. . ,m(n) with 

f T f, lA„, fc (t,£,a)dtd£K^)<oo, 

such that G n converges to G in L 2 ([0,T] x O x R, dsd£v(da); H) if n tends to infinity. Note 
that the processes 

I / / G(s, £, a) Tt(ds, d£ } da) 1 and ( / / G n (s,!:,a)7r(ds,dt:,da) 1 , 

VJo JOxR /te[0,T] VJO JOxR /t€[0,T] 

n e IN, are martingales and hence have cadlag modifications. Applying Doob's inequality for 
submartingales, we have 

I rt r rt r 2 

Ee 



sup / / G(s,£,a)-rt(ds,df,da)- / / G n (s,£,a)Tr(ds,d£,dcr) 

t< T \J0 JOxR Jo JOxR 

<4E| f T f (G(s,Z,<T)-G n (s,Z,cT)) TT(ds,dC,da) 
\J0 JOxR 
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H 



and the last term goes to zero as n tends to infinity. 

Remember that the Poisson random measure tt on [0,T] x O x R with given compensator 
dtd^v(da) is given by 



7T= £%^,E 3 ), 
.7 = 1 



where {(Tj, Sy, Xj))_. £l>J is a properly chosen sequence of random elements in [0, T] x x R, cf. 
|27j . Chapter 6. Since the functions G n , re e IN, are simple functions, the stochastic integral of 
G n with respect to n is just an cu-wise integral with respect to 



E 6 (Tj (w),^ (w)) (rf* 3 d£,da)-dt d£ v(d<j) . 



Therefore, for every re e IN", the cadlag process 

( r f G n ( S ^,a)7r(d S! dC,^)) 



te[0,T] 



has the jumps G n [Tj(oo),'Ej(ui),'Ej(uj)^, j e IN, Tj(ijj) < T, occurring at the jump times 
Tj(uj), j e IN, Tj(oj) < T, for almost every uj e Q. 

The convergence G n -*■ G, n -» cxj, in L 2 ([0, T] x0xR, dtd^v{da)] H) implies the convergence 
G nj (t,^,a) -> G(t,^,a), j -*■ oo, in iT for dsd£i/(dcr)-almost every (t,£,cr) e [0,T] x x R 
(where (G n . )j £ fj is a subsequence of (G n )nelN)- This, the uniform convergence following from 
(lA.ip and the fact that the laws of the random vectors (Tj,Sj,'Sj), j e M, are absolutely 
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continuous with respect to dtd^v{d(j) imply that, for almost every uj e £1, the jumps of 
the process {Yh,t) t€ rQ T -i occur at the jump times Tj(uj), j e IN, Tj(ui) < T and are exactly 

G(T>), j e IN, Tj(u) < T. 
Consequently, we have 



E 



! Y, {MT - 1, Y h , t ) - v h (T - 1, - (D x v h (T - 1, Yh,t-), AY^) } 

t<T 

= E £ {^{T-Tj^T^ + GiTj^Ej^-vniT-Tj^T..) 
jeN,Tj<T 

- {D x v h (T - Tj,Y htTj -), G(T j ,E j ,'E j )) H }. 
With similar arguments, one can show that 

E f T f {v h (T-t,Y h j_ + G(t,£,a))-v h (T-t,Y hit _) 

Jo JOxK I v ' ' 

- (D x v h (T - t, Y h}t _), G(t, e, a)) R } n(dt, d£, da) 

= E Y {vhiT-T^Y^ + GiTj^j^-VhiT-Tj^.) 

- {D x v h (T - Tj,Y h ^), G(Tj,Ej,Y^j)) H }. 

This finishes the proof of ( 14 0|) , 



B Extension to the case /{ui>i} cr 2 v(da) = oo 

Here we show how to handle the case where assumption © fails to hold, assuming just that 
K{0}) = and 

f mm(a 2 ,l)v(da) < oo. 0) 

(Zf)te[o,T] as defined in Section 2 is a process with values in a Hilbert space U, which we 

assume to be of the form U - H~z~ e (0), e > 0, in this section. Informally, Z% is also denoted 
by 

Z(t,d£) = / / a Tr(ds,d£,da). 

As we cannot follow the construction in Section 2 if assumption ([7J does not hold, we consider 
in this section a process of the form 

Z(t,d£)= / / <jir(ds,d£,d<j) + / / a ir(ds,d£,da) 
Jo J{M<i} Jo J{\a\>l} 

=:M(t,d£) + P(t,d€). 

More rigorously, we define the process {Zt)te[o,T] by Zt '■= Mt + Pt, t e [0,T], where (Mt)tn[o,T] 
is an impulsive cylindrical process with jump size intensity measure UM(dcr) = \[_i i](o)v{d<j) 
and (Pt)tt[o,T] i s a ^-valued compound Poisson process given by 

Tj<t 
Pil>l 
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Here ((7^,Ej-,Ej)) . g]N is the sequence of random elements in [0, T] x O x R introduced in 
Section 2 and for (eOwe denote by 5^ e U the Dirac delta function at £. 
The equation 

dXj + = Q 1/2 (iZt = Q 1/2 (dM t + dP t ), X = x £H, te[Q,T], gf) 

has a weak solution with values in H if the linear operator Q 1 ^ 2 ■ H -*■ H has a (not necessarily 
bounded) extension Q 1 ! 2 : U -* H, compare Example [10] below. In this case a weak solution 
to (jlT]) can be constructed as follows (compare |27] . Chapter 9.7). Let us consider a sequence 
of stopping times 

r m := inf {t 6 [0, T] : P t - P t _ f V m }, m e M, 
where the Vm are bounded measurable subsets of U given by 

V m ■= {oS^ : £ e O, a e [-m,m]} , m e IN, 

and inf0 := oo. Note that due to (ITfl) we have P(Umew{ r m = °°}) = 1- For every m e M we 
define 

P^W^^), M t ™:=P™-iu m , t 6 [0,T], 
with u m := EP| m -' e J7, and consider the problem 

dx\ m] + Ax\ m] dt = Q 1/2 ( u m dt + dM t + dM\ m] ) , X^ m] = x e ff. fll') 

(M ( ^ is an impulsive cylindrical process in the sense of Section 2, with jump size 

intensity measure l[- TO) m]\[-l,l]( cr ) v(da\ There exists a unique weak solution (jfj ) t6 r yi 
to (HQ) for every meK, given by 

X t [m] = S(t)x + ^ S(t - s)<5 1/2 («m ds + dM s + dMl m] J, t e [0, T]. 

For m < n e M one has {r m = oo} c {r„ = cxd}, and for every t e [0,T] the equality 
Xj m] l {rm=00 } = Xj n] l {rm=00 } holds P-almost surely. Consequently, the process (^t) te [ ,T] 
given by the P-almost surely unique limit 

X t := lim 4 m] l {Tm=oo} 
is a weak solution to equation (HTjl . 

Example 10. Equation (Hfj) has an i?- valued weak solution (Xt)te[o,T] if t ne operator Q 1 ! 2 
can be extended to Q 1/2 :U^H. Here we give an example for such an extension. Consider 
an operator Q as in the example mentioned in the introduction, i.e. 

Q 1/2 x(0 = f Q <7o(£, CMC) dC, xzL 2 (0),ZeO, 

where qo is a positive semidefinite symmetric function on 0x0. Assume furthermore that 
Q l l 2 e L(H, U*) = L(L 2 (0),HQ /2+e (0)). Then Q 1 / 2 can be extended to a linear and bounded 
operator e L(U,H). For u<iU and x ^ H set 

{Q 1/2 u,x) H := v ,(u, f qo(- ,v) x (v)dv) v = u *(u,Q 1/2 x) u , 
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v * (■ , - )u '• U* x U -*■ R being the canonical dual pairing. By the Riesz representation theorem, 
this defines Q 1 ! 2 - U -*■ H. In particular, for £ e O we have the equality in H 

Q V % = qo(-,0- 

Remark 11. If the extension of Q 1 ! 2 is continuous, i.e. if Q^^ 2 € if), then the term Qu m 
appearing in equation (}4Tj) can be represented as follows: Let /i m denote the jump intensity 



measure of (pj- m ^ ) t6 r Q T -i ■ We have 

u m = z[j, m (dz)= / aScd^u(da) 

JU JOx([-m,m]\[-l,l]) 

JOx([-m,m]\[-l,l]) [ fe=1 J 

where all integrals are [/-valued Z^-Bochner integrals, (ek)ke¥l is an orthonormal basis in 
H and the infinite sum is a limit in the space L 2 ((D x ([— m,m] \ [-1, 1]), d£ u(da) ; t/J c 
L\0 x ([-m,m] \ [-l,l]),d£z/(do-); U). 
Consequently, 



fc=1 L->Ox([-m,m]\[-l,l]J J 

the infinite sum being a limit in H. 

We extend the approximation scheme (|20p - (|22p of Section 4 and define discretizations 
(X^' )ne{o,...,JV} °f (^t)t€[o,T] f° r integers N, m > 1 and real /i > by 

n-1 

fc=0 

where n e {0, . . . , N}, compare (|24p . Given a bounded function 99 : H -> R we find 

< \E<p(X^' [m] ) -E^(4 m] )| + |E^(4 m] ) -Ev9(X T ) 

< |E</?(xf ,[m] ) - E^(x[ m] )| + 2|M| 00 P(r m < 00) 

< |E</>(xf ,[m] ) - E^X^)] + 2HIOO min(rA(C) i/(R \ [-m,m]) , l), 

where ||^|oo : = su PxeH I^C^OI-ff an d ^(C) i s the Lebesgue measure of O. Of course, due to 
(fTfj) . z^(R s [— m, m]) tends to zero as m -*■ 00 and the order of convergence depends on the 
particular jump size intensity measure v. 

The term |E^(xf '^-E^X^)] can be treated with the methods developed in Sections 
4, 5 and 6. We briefly sketch the alterations. Let 



X]™ ] = S h (t)P h x + J Q S h (t-s)P h Q 1/2 (u m ds + dM s + dMl m l), t€[0,T], ^M) 
Y lf S h (T)P h x + (u m ds + dM s + dM™) , t 6 [0,T], ([26l) 



25 



S% At P h x + £ T( S ) (u m ds + dM s + dM^ ) , f e [0, T], (EH) 

with = S h (T - s)P h Q 1 / 2 and r(s) = E^^^r^Atl^^jC^J^Q 1 ^ as in Section 

4. Let u[ m] : [0,T] x F ^ R and : [0,T] x if R be defined by 

(t, x) := E<p (x + Jj t S h (T - r)P h Q x < 2 (u m ds + dM s + dM™ ) V (EI) 
v [m] (t,x) := Eip^x + Jj ^(T - r)Q 1/2 (u m ds + dM s + dMl" 1 ^. QM) 
Arguing along the lines of Section 5, we get 

Theorem EJ. Assume that (J2Q) holds, \\S(t)Q 1/2 \\( HS) € L 2 ([0,T],dt) and that Q 1/2 has a 
linear and continuous extension Q 1 / 2 : U -+ H . Let ip e C 2 (H) such that for each x e H the 
derivative D 2 ip(x) is an element of L( HS j(H) and that the mapping x D 2 ip(x) € L( H gj(H) 

is uniformly continuous on any bounded subset of H. Let T > 1, m € IN and (^Q )t € [o,T] 
be the H -valued stochastic process satisfying equation §4\'\i . For any N > 1 and h > 0, let 

(Xfo' )ne{o,...,Ar} ^ e ffiwen fry (|24ip and iei (-Xjj^ )te[o,T] be as in f)25f ]) . 
Then, the following error expansions hold for the time discretization error 

E^X^-E^X^) = {vt ] (T,S» At P h x )-vl™\T,S h (T)P h xo)} 

+ (a) ( A4 m] (T - t, ) , a) - G(t, M) H \dt <% u{da) , 

and for the space discretization error 

E^(4™ ] )-Ev9(4 m] ) = {v^(T,S h (T)P h x )-v^(T,S(T)x )} 

JO JOx[-m,mJ x ' x ' 

+ 1 [_!,!] (a) ( (T - i, ) , £(t, e, (r) - F(t, M)\dt d£ v(do-). 

Remark 12. In some cases the assumption ||5 , (t)Q 1 / 2 ||(HS) e L 2 ([0,T],dt) in Theorem [Sf 
is already fulfilled as a consequence of Q 1 ! 2 having a linear and continuous extension Q 1 ! 2 ■ 
U -> H, because the latter is equivalent to Q 1 ! 2 € L(H,U*) = L(L 2 (O) , H^ 2+t (0)) . As an 
example, let O = (0, l) d and (A,D(A)) = ( - A, #£(0) n H 2 (0)). Setting /3 := § + e, we 
have A^^Q 1 ! 2 e L(iJ) due to Q 1 / 2 e L(H,U*) and this yields |£(*)Q 1/2 |(HS) e ^ 2 ([0,T],dt) 
because Tr(74~^) < oo. 

If the jump size intensity measure v satisfies 

min(|cr|, 1) v{da) < oo, (H) 
one gets the following counterpart of Theorem [8] 



JR 
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Theorem M • Assume ® and (JTDJ) ? JTTJI, d), CHI), CGI listed on page [731 Assume /ur- 
thermore that Q l l 2 has a linear and continuous extension Q l l 2 ■ U -> il. Let 95 € C 2 (H) such 
that for each x € H the derivative D 2 ip(x) is an element of Lfugj(H) and that the mapping 
x >-> D 2 ip(x) e L( HS j(H) is uniformly continuous on every bounded subset of H. Let T > 1 and 
(-X^)t 6 [o,T] &e the H-valued stochastic process satisfying the equation 

dX t + AX t = Q 1/2 dZ 4 = Q 1/2 (dM t + dP t ), X = x €H, te [0, T]. 

For any TV", m > 1 and /i > 0, let (X^ m ^) ne { jj\ ^ e 5™ eri % ([241 p . 

T/ien /or a//7<l-a + /3<l t/iere exists a constant C = C(T, <p, A, Q, \xq\, 7) > which does 
not depend on h, N and m such that the following inequality holds 

\M^ [m] )-^{x T )\ 

\a\ v(do)(h 21 + (At) 7 ) + 2|^||oo min (TA(O) i/(R \ [-m,m]) , l), 

m 

w/iere At = T/N < 1. 
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